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Abstract. We suggest that non-thermal emission can be produced by multiple scatterings of the photons between the 
supernova ejecta and pre-shock material in supernova shock breakout. Such bulk-Comptonization process may significantly 
change the original thermal photon spectrum, forming a power-law non-thermal component at higher energies. We then show 
that the luminous X-ray outburst XRO081009 associated with SN2008D is likely to be such shock breakout emission from an 
ordinary type Ib/c supernova. 
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INTRODUCTION 

Supernova shock breakout has been predicted for a few decades (e.g. |H]; JUl; 01; J3l). In core collapse supernovae, a 
shock wave is generated which propagates through the progenitor star and ejects the envelope. As the shock propagates 
through the envelope, it is mediated by radiation: the post shock energy density is dominated by radiation, and the 
shock transition is mediated by Compton scattering. As the shock approaches the edge of the star, the optical depth 
of the plasma lying ahead of the shock decreases. At the radius where the optical depth drops to ~ c/v s , where v s 
is the shock velocity, the shock wave emerges through the surface and the radiation decouples, accompanied by a 
very bright ultraviolet/X-ray burst of radiation. The radiation spectrum is expected to be thermal-dominated. The term 
"shock breakout" is commonly used to refer to the emergence of the shock from the edge of the star. However, if the 
star is surrounded by an optically thick wind, the radiation mediated shock would continue to propagate into the wind, 
up to the point where the wind optical depth drops below ~ c/v s . In this case, shock breakout, occurs as the shock 
propagates through the wind, at a radius which may be significantly larger than the star's radius. We use here the term 
"shock breakout" to denote this transition in general, regardless of whether it occurs as the shock reaches the edge of 
the star or further out within an optically thick wind. 

This kind of radiation has previously never been directly detected from any normal supernovae due to its transient 
nature and its very early (minutes to hours) occurrence, in the absence of a suitably prompt trigger alert. On 2006 
February 18, thanks to its sensitive gamma-ray trigger and rapid slewing capability, Swift has detected early thermal 
X-rays emission from a supernova (SN2006aj) associated with a low-luminosity GRB, namely GRB060218 (e.g. J5tl), 
which has been interpreted as arising from the breakout of a radiation-dominated shock (@1; (^). The long duration 
(~ 3000 s) of the thermal X-ray emission of GRB060218/SN2006aj suggests that the shock breaks out from a dense, 
optically-thick wind surrounding the progenitor staiQ. The energetic of the thermal and non-thermal emission is of the 
order of 10 49 ergs, which is much larger than the predicted radiation energy from a normal type Ib/c supernova (|01). 
Such a large energy release could be due to a larger kinetic energy in hypernova SN2006aj and/or an central engine 
that drives a jet which is, however, chocked in the outward propagation. 

On 2008 January 9, a bright X-ray outburst XRO080109 was serendipitously discovered during a scheduled Swift 
observations of the galaxy NGC 2770 ([8]). An ordinary type Ib/c supernova in coincident with this outburst was later 
spectroscopically identified and named SN2008D. This outburst has an energy of 2 x 10 46 erg, which is three order 



Another possibility is that shock breaks out from an optically thick shell pre-ejected from the progenitor (the required shell mass is only 1(T 7 M , 
seed). 



of magnitudes smaller than even low-luminosity GRBs, but is astonishingly close to the predicted shock breakout 
radiation energy from a normal type Ib/c supernova (01). 

The most mysterious thing of this X-ray outburst is the non-thermal spectrum, which is in contrast with the thermal 
spectrum predicted by the shock break theory. There is a general agreement that a power-la w sp ectrum (with photon 
index T — 2.3 ± 0.3) provides a better fit for this x-ray outburst than a blackbody (Hfl; [Ml; lIlOll : lUllO . The temporal 
evolution is characterized by a fast rise and exponential decay, with a FWHM duration about 100 s. The rise time is 
fitted to be 63 ±7 s (H). 

NON-THERMAL EMISSION FROM SUPERNOVA SHOCK BREAKOUT DUE TO 

BULK MOTION COMPTONIZATION 

Blandford & Payne [12] first noted the importance of bulk motion acceleration of photons in a radiation-dominated 
shock. They found that photons are preferentially upscattered by the bulk motion rather than by the thermal motions of 
the electrons, and a power-law spectrum extending to high energies forms, when the electron thermal velocity is less 
than the shock velocity v s . Repeated scatterings using the energy of the bulk motions of two approaching relativistic 
shells in the context of GRB internal shocks was studied by [ 13]. They found that the seed synchrotron photons can be 
boosted to much higher energies, which is confirmed by their Monte Carlo simulations. This process is equivalent to 
the Fermi acceleration mechanism of particles or photon scattering off Alvfen waves ( ifTill ). but here the mechanism, 
instead, uses the relative bulk motion and accelerates photons. 

Let's consider a mildly relativistic ejecta driving a radiation-dominated shock into the stellar envelope of supernova 
progenitor or into an optically thick wind (or pre-ejected shell) surrounding it. This mildly relativistic ejecta could 
result from the shock acceleration in the surface of the type Ib/c supernova progenitor ([4]) or from a jet that is 
chocked and expands sideways in the outward propagation inside the star. Once the optical depth of the material in 
front of the shock drops below c/v s (where v s is the shock velocity), the photons escape and produce a breakout flash. 
Since the Thompson scattering optical depth is non-negligible in front of the shock while the shock is breaking out, 
some fraction of the thermal photons will be scattered back. The back-scattered photons will be scattered forward by 
the expanding ejecta or shocked plasma, boosting up their energy. The backward-forward scattering cycle may repeat 
itself many times for some fraction of the photons, boosting their energy by a large factor J3]. 



A qualitative description 

Ignoring the time dependence of t, the optical depth ahead of the shock, the physical situation is similar to that of 
Comptonization by a thermal electron plasma. In many cases, the electrons are cold and their momentum is dominated 
by the bulk motion. Assuming each scattering amplifies the photon energy by a factor A, the energy of a photon 
escaping after k scatterings is = £jA , where e, and are the initial and final photon energies respectively. For fixed 
T, a photon scattered by the ejecta has a probability 1 — e~ T to be scattered back towards the ejecta, and a probability 
e~ T to escape. The probability for a photon to undergo k scatterings before escaping is (1 — e~ z ) k , and since the photon 
energy is multiplied by A per scattering, the escaping photon intensity would have a power-law shape 

F(e k ) ~Ffe)(l -0*~^i)(£*Mr a (1) 

with 

a = -ln(l-e- T )/lnA (2) 

The photon energy amplification factor A is determined by the kinetic energy of the electrons. For trans-relativistic 
electrons and isotropic photon distributions, A ~ T 2 (l + j3 2 /3). The power-law spectrum extends to a cutoff energy, 
which is the smaller of the electron kinetic energy, ~ (r — 1 )m e c 2 , and its rest mass m e c 2 (due to the Klein-Nishina 
effect). 

An important difference between the usual thermal electron (or bulk) Comptonization case and the current case is 
that in the present case the scattering optical depth decreases with time as the mildly relativistic ejecta moves outward. 
Initially, when the optical depth T > 1, the slope of vF v is positive and most radiation is emitted at high energies. As T 
decreases, the spectrum becomes softer and softer, and at late times the spectrum is composed of a thermal peak plus 
a weak high-energy power-law tail. In general, we expect a noticeable spectral softening of the nonthermal emission 



-10 12 3 

Log(hv/KeV) 



FIGURE 1. The time-integrated energy distribution of the escaping photons. 10 6 photons, with a black body distribution at 
kBT/h = 0. 15KeV (black dotted line), are injected at four different times, corresponding to optical depths of the pre-shocked medium 
of Tj n j = 2, c/v s , 0.5 and 0.2. The ejecta Lorentz factors are T = 1.5 and T = 2 for upper panel and lower panel respectively. Note, 
that the "humps" seen in the spectra are an artifact of the one-dimensional simulation, and are expected to be smoothed out in 
reality. See |7] for more details. 



with time. This spectral softening is expected to be accompanied by a decrease in the Compton luminosity. At early 
time, when the effective Compton parameter Y =A(1 — e~ T ) > 1, the Compton luminosity may exceed the thermal 
luminosity (it is limited by the kinetic energy of the ejecta E^). We expect the Compton luminosity to decrease with 
time, as Y decreases. 

The x-ray or gamma-ray light curves produced in this model generally have a simple profile without multi-peak 
structure. The characteristic variability timescale 8t of the burst is determined by the the radius R where the optical 
depth of the material ahead of the shock drops to ~ 1, i.e. 8t ~ R(t = l)/c (If the stellar wind surrounding the 
progenitor were optically thin everywhere, the shock would break out from the SN progenitor stellar envelope and the 
variability time would be about R*/c, where is the stellar radius.) 



Monte Carlo simulation of photon "acceleration" 

In order to understand the photon "acceleration" mechanism in the time-dependent case, we carried out a Monte 
Carlo simulation of repeated Compton scattering during shock breakout from a dense stellar wind [7] . We approximate 
the hydrodynamics of the problem as follows. We consider the mildly relativistic ejected shell to act as a piston with 



a time-independent bulk Lorentz factor F and an infinite optical depth T e j. Since the stellar wind swept up by the 
ejecta is not sufficient to decelerate it, the constant velocity of the interface is justified. This"piston" drives a shock 
into the surrounding medium, where the density profile is assumed to follow n °< R- 2 . We assume the shock width 
to be infinitesimal. There are three distinct regions in this picture: The moving piston, the shocked medium and the 
pre-shocked medium. The shocked medium is considered to form a homogeneous shell, and the electrons are regarded 
as cold, with a bulk velocity same as that of the ejecta. The velocity of the shock front can then be obtained consistently 
from the shock jump condition. We further simplify the problem by considering a one-dimensional situation, where the 
motion of the photons is confined to one dimension, forward or backward relative to the shock expansion direction. 
The photons are injected at the shock front and would then be repeatedly scattered among these three components 
until they escape out to a sufficiently further region. Our simulation takes into account the proper possibility of photon 
scattering by each of the three components (with the appropriate photon energy gain or loss). The photon scattering 
probability by the unshocked wind and shocked wind is determined by the optical depth of each component, while the 
piston is regarded as a mirror since it has a larger optical depth much larger than unity. 

We study the spectra of escaping photons, which result from the "injection" of photons at the shock front at various 
radii, corresponding to various values of T at the injection time, denoted by T; n j. The photons are treated as "test 
particles", and their scattering history is followed as the shock expands and T decreases. The resulting time-integrated 
vF v spectra are given in Fig. 1 for two values of T, F = 2 and F = 1.5. The red, blue, olive and black curves describe 
the spectrum of escaping photons resulting from the injection of a thermal distribution of photons at T; n j = 2, c/v s , 0.5 
and 0.2 respectively. The injected photons are assumed to have a thermal spectrum with T = 0.15 keV. The "humps" 
seen in the simulated spectra correspond to different orders of Compton scattering, with two nearby humps separated 
by ~ r 2 . These humps appear because we do not consider in our one-dimensional simulation the angular distribution 
of scattered photons. In reality, these humps are expected to be smoothed out. Fig. 1 demonstrates that a significant 
fraction the thermal photons that are injected at optical depth T; n j ~ 1, where the photons are expected to escape the 
shock, may be "accelerated" to high energy. The resulting non-thermal component may carry a significant fraction 
of the shell energy, and its luminosity could exceed that of the thermal component, depending on the value of F. As 
expected, the spectrum depends on F and on Tj n j , with harder spectra obtained for larger values of F and T; n j . Note that 
the spectra are steeper (softer) at high energies, due to the longer time required for acceleration to higher energies, 
which implies a significant decrease in t during the acceleration process. 

SHOCK BREAKOUT ORIGIN OF THE X-RAY OUTBURST XRO 080109 

Although the thermal x-ray component is not identified in XRO 080109, there are a few facts that favor the shock 
breakout interpretation: 1) smooth light curve-consistent with shock breakout prediction; 2) low total radiation energy 
of 10 46 erg-close to the theory prediction for a normal type Ib/c SN; 3) early (few days) decline in the optical emission 
from SN, probably optical emission from the adiabatically cooling ejecta heated by shock breakout (la]). 

The non-thermal spectrum of XRO 080109 can be due to the bulk-comptonization of the thermal shock breakout 
photons. Then, the peak of the thermal component should be not larger than 0.3Ke V (the low threshold energy of XRT) 
so that the thermal component was not detected by XRT. Since the Compton Y parameter is typically of the order of 
unity, the energy in non-thermal emission should be roughly comparable to that in thermal energy. From these simple 
arguments, we now derive the constraint on the shock parameters, following @]. 

For a strong radiation-dominated shock, the radiation pressure p = f(Ffi) 2 pc 2 , where p is the mass density of the 
pre-shocked density and / ~ 0.8 for trans-relativistic shock (la]). The post shock temperature T c [ is related to the the 
postshock pressure by aT^ = 3p, and the observed temperature is T = JdTj, where Yd is the Lorentz factor of the 
downstream flow. The temperature of the shock breakout thermal emission T is related to the velocity of the shock 
and photosphere radius by 

ar* = 3/(r/3)VM - (r/3) V 2 = (r/3) 6 c 2 (-U, (3) 

KK ph 

where the optical depth t(R) = KRp with T(Ri, r ) = t(R p i,) = 1 at the breakout radius Rt, r for trans-relativistic shock 
and K is the Thomson opacity. 



Assuming the energy in the thermal component is comparable to that in the non-thermal one (which is likely for a 
trans-relativistic shock), i.e. E t h ~2 x 10 46 erg, we can obtain, from Eq.(2) in ref.fl, 
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and 



(TBY — = (4) 

1 P) "0.5x4^ 2nR 2 ph c 2 K ' 

Combining Eqs.(3) and (4), we obtain 

TR = 1 l( T ) 2/7 ( ^ ) l ' 2 \ - ) 3 / 28 (5) 

P [ 0AKeV } l 2xl0 46 erg J { 0.2g- 1 cm 2) K) 

"" = i* 10,, ^>- 4/7 <27^» 3/, (a2^) 2, '- < 6 > 

where K ~ 0.2g _1 cm 2 is the Thompson opacity for ionized He wind. 

Requiring T < 0. IKeV to account for the non-detection of the thermal peak, we derive the constraints on the shock 
velocity Tj3 and breakout radius R/, r , i.e. 

TB < 1 U — )'/ 28 ( - ) 3 / 28 (7) 

P " l 2xl0 46 erg J [ 0.2 g -^cm 2> ' K) 

R br > 1 x 10 n ( . f* ) 3/7 ( * 2 ) 2/? cm. (8) 
2 x 10 4O erg 0.2g 'cm z 

We note that there is a misuse of the relation L rac i = 4nR 2 aT 4 to derive the shock breakout radius in some papers, 
where T is the radiation temperature. This is because the light travel time ( R/c) will lengthen the breakout emission 
duration significantly and render the observed radiation luminosity lower than 4%R 2 aT 4 significantly (see also |01). 
So one would get a much smaller radius for shock breakout if this incorrect relation is used. 
The rise time of the x-ray transient is about 60 s ([8], implying a radius for the emitting region 

R ph ~ c8t = 2 x 10 12 cm, (9) 

which is consistent with the above inferred breakout radius. Since this radius is larger than that of the WR progenitor, 
the shock must break out from the surrounding stellar wind. From T = 1 at Rt, r , we derive M — 4nv w Rt, r / K > 
7 x lO^Mcgyr-^K/O^g^cm 2 )- 1 for v w = lOOOKms- 1 . 

So we need a trans-relativistic shock with TB < 1 at the time when shock is breaking out. This can be achieved for 
a normal type Ib/s SN explosion, according to the calculation of fl. Substituting the inferred ejecta kinetic energy 
(E K = 2 -4 x 10 51 erg) and mass (M ej = 3 -5M ) the of SN 2008D and a typical Wolf-Rayet star radius ofR = 10 u cm 
into Eq.(32) of [4], the maximum velocity that the radiation-dominated shock can reach is trans-relativistic. 

With a velocity 0.5 < T/3 < 1 , the energy amplification factor for one Compton scattering is 

A = r 2 (l+j3 2 /3) = 1.5-2.3 (10) 
which implies an equivalent Compton parameter Y = A(l — e~ T ) ~ 1, consistent with our earlier assumption. 



DISCUSSIONS 

Chevalier & Fransson fl5ll questioned the bulk Comptonization mechanism for XRO080109/SN2008D by arguing that 
"the breakout radiation is capable of accelerating the matter ahead of the shock front so that the formation of the gas 
dominated shock is delayed". This is not true because as long as there is a converging flow (in our case the supernova 
ejecta moving relative to the pre-shock matter), the repeated scatterings will occur and the bulk Comptonization will 
work. So the formation of an viscous shock is not a necessity for the bulk Comptonization process. Also, for a type 
Ib/c supernova shock that reach a mildly-relativistic velocity, the radiation is actually incapable of accelerating the 
matter ahead of the shock front to such a mildly-relativistic velocity, according to our calculation given below. The 
radiative acceleration is 

5 4tzR 2 c y Q.2s- l cm 2jy lQ 42 eTg&- 1 ' 12 V ; 



So the maximum velocity that the matter ahead of the shock front be accelerated to is 

v m = gR At = 5 x 10 9 ( \ , )( .J r " d T )R^(-^-)cm S -\ (12) 
0.2g l cm z 10 4Z ergs 1 lL 100s 

which is below the mildly-relativistic velocity of the shock front. 
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